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The top quark is the heaviest known fundamental par- 
ticle in the standard model of particle physics. Since the 
1995 discovery of the top quark at the Fermilab Teva- 
tron [T], both the CDF and DO experiments have been 
improving the measurement of its mass m t , which is a 
fundamental parameter in the standard model [2]. Loop 
corrections in electroweak theory relate m t (along with 
the W boson mass m\y) to the mass of the predicted 
Higgs boson. Thus, precision measurements of mt help 
to constrain the value of the Higgs boson mass [3]. 

This Letter describes the single most-precise measure- 
ment to date of the top quark mass. It is performed 
on data collected by the CDF II detector [4 during 
Run II of the Fermilab Tevatron pp collider operating 
at y/s — 1.96 TeV with a total integrated luminosity of 
5.6 fb _1 . The measurement is performed on candidate 
ti events containing a lepton and four jets J5L For each 
event selected in this analysis, we calculate the probabil- 
ity of observing that event by integrating the matrix ele- 
ment for ti production and decay over phase-space vari- 
ables. We use a neural network discriminant to distin- 
guish between signal and background events to correct 
for the contribution due to background, and employ a 
cut on the peak likelihood for a given event for addi- 
tional rejection of background and poorly modeled sig- 
nal events. This analysis offers a gain of nearly 20% in 
statistical precision over our previous measurement [3] 
given an equal number of events; there, in order to make 
the likelihood integration computationally tractable, we 
introduced kinematic assumptions to reduce the dimen- 
sionality of the integral. In the present analysis, we use a 
quasi-Monte Carlo integration technique [6_, which con- 
verges more rapidly than the typical 0(N~ 1 ^ 2 ) conver- 
gence of standard Monte Carlo integration. This allows 
us to integrate over a total of 19 dimensions in a com- 
putationally practical time, resulting in a more accurate 
modeling of the event. Furthermore, in addition to the 
increased data sample available with more integrated lu- 
minosity delivered by the Tevatron, we have expanded 
our muon selection ability, which increases the size of 
our data sample by nearly 30%. In total, this measure- 
ment improves our statistical precision by a factor of two 
over our previous analysis with 1.9 fb _1 [5]. 

In this measurement, the largest uncertainty is due to 
the uncertainty in the jet energy scale (JES) determina- 
tion. To reduce this uncertainty, we calculate the likeli- 
hood as a two-dimensional function of mt and a second 
parameter, Ajes, which corrects the jet energies by a 
factor of 1 + Ajes 'Ojj where o.j is the fractional system- 
atic uncertainty on the energy for a given jet 0[8]. The 
known W boson mass is used to constrain the W — > qq' 
decay, which yields information on the Ajes parameter. 
We can thus optimally combine events to reduce the total 
uncertainty on m t due to JES. 

Within the standard model, the top quark decays al- 



most exclusively into a W boson and a b quark. We 
define a "lepton -I- jets" event as an event where one 
of the W bosons produced by the ti pair decays into a 
charged lepton (in this analysis, an electron or a muon) 
and a neutrino, and the other into a qq' pair. The two b 
quarks and two quarks from the W boson then produce 
jets in the detector [§]. We thus require candidate events 
to have an electron with Et > 20 GeV or a muon with 
Pt > 20 GeV/c in the central detector (\rj\ < 1 [TO]), or a 
muon with pt > 20 GeV/c obtained with a trigger [3] on 
missing transverse energy, Et [ID] , instead of a central 
muon. As the neutrino energy is not detected, we require 
Et > 20 GeV in the event. We also require exactly four 
jets with Et > 20 GeV within the region \rj\ < 2.0, at 
least one of which must be tagged as a & jet using a sec- 
ondary vertex tagging algorithm |llj . To model ti events, 
we use Monte Carlo-simulated events generated with the 
pythia [12] generator for 15 different m t values ranging 
from 162 to 184 GeV/c 2 . 

Background events contributing to the selected sam- 
ple are: a) events in which a W boson is produced in 
conjunction with heavy-flavor quarks (bb, cc, or c); b) 
events in which a W boson is produced along with light 
quarks, at least one of which is mistagged as heavy fla- 
vor; c) QCD events that do not contain a true W bo- 
son; d) diboson (WW, WZ, or ZZ) or Z + jets events; 
e) single top events. We model the contribution from 
W+jets events using ALPGEN [T3], single top events using 
MADGRAPH [2], and diboson events with pythia. The 
Z+jets contributions are not modeled separately, but are 
included in the W+light flavor contribution. All Monte 
Carlo samples are processed with the CDF II detector re- 
sponse simulation package [TS]. The non-IV QCD back- 
ground is modeled using a sideband of data events se- 
lected to have a small contribution from heavy boson 
decay. The numbers of background events are estimated 
with the method used for the tt cross section measure- 
ment [16] . and are shown in Table [i] 

For each event, we construct a likelihood as a function 
of m t and Ajes using the following integral: 



L(y\ m u Ajes) 



1 



1 



N(m t )A(m t , Aj ES ) 

24 

^2wiLi(y\ m t , Ajes) 
i 

? 

\M{m u x)\ 2 d$(5), 



(1) 



L%{y | m t , Ajes) = / - y ~ F p~" TF(/7 VA. ILS )x 



where y are the quantities measured in the detector 
(the momenta of the jets and charged lepton), x are 
the parton-level quantities that define the kinematics 
of the event, N(m t ) is a global normalization factor, 
A(mt, Ajes) is the event acceptance as a function of 
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TABLE I: Expected sample composition for an integrated lu- 
minosity of 5.6 fb _1 . The ti contribution is estimated using 
a cross-section of 7.4 pb [17] and mt = 172.5 GeV/c 2 . 



Event type 


1 b tag 


> 2 b tags 


W+heavy flavor 


129.5 ± 42.1 


15.7 ± 5.5 


non-W QCD 


50.1 ± 25.5 


5.5 ± 3.8 


IT+light flavor mistag 


48.5 ± 17.1 


1.0 ± 0.4 


diboson (WW, WZ, ZZ) 


10.5 ± 1.1 


1.0 ± 0.1 


Single top 


13.3 ± 0.9 


4.0 ± 0.4 


Z ™> U + jets 


9.9 ± 1.2 


0.8 ± 0.1 


Total background 


261.8 ± 60.6 


28.0 ± 9.6 


ti signal 


767.3 ± 97.2 


276.5 ± 43.0 


Total expected 


1029.1 ± 114.5 


304.5 ± 44.1 


Events observed 


1016 


247 



m t and Ajes, /(zi) and f(z2) are the parton distri- 
bution functions (PDFs) for incoming parton momen- 
tum fractions z\ and z%, FF is the relativistic flux fac- 
tor, TF(?7 | x, Ajes) are the transfer functions that de- 
scribe the measured jet-momentum distributions given 
the quark kinematics, the phase space for the 

eight particles in the ti production and decay process, 
and M(mt, x) is the matrix element for the process. The 
integral is calculated for each of the 24 possible permu- 
tations of jet-parton assignment and then summed with 
weights Wi determined by the probability that a b or light 
parton will result in a ^-tagged or untagged jet. 

We use the Kleiss-Stirling matrix element [18], which 
is a leading-order matrix element including both qq — > ti 
and gg — > ti production processes, as well as all spin 
correlations. For the PDFs, we use the CTEQ5L func- 
tions |19j for the incoming qq and gluons. The normaliza- 
tion factor N(m t ) is obtained by integrating the Kleiss- 
Stirling matrix element with the PDFs and the flux factor 
over the phase space formed by the two initial and the 
six final-state particles. The acceptance A(m t ,AjEs) is 
obtained from simulated events where the parton direc- 
tions and momenta are smeared to simulate final-state 
jets. The transfer functions connect the measured jets to 
the partons. We construct the transfer functions by tak- 
ing simulated ti —> lepton + jets events in a wide range 
of masses and matching the simulated jets to their parent 
partons. The transfer functions are separated into mo- 
mentum and angular terms; both are constructed with 
dependence on the true jet pt and mass from the Monte 
Carlo simulation. The transfer functions are constructed 
separately for b and light quarks, as well as for each of 
four bins of jet rj. There are 32 phase space integration 
variables in Eq. ([I]) (for the two initial partons and six 
final partons). Four of these are eliminated by energy 
and momentum conservation, and four more by taking 
the charged lepton, neutrino, and initial parton masses 



as known. In addition, we assume that the lepton mo- 
mentum is perfectly measured, and we neglect the effects 
of the individual transverse momenta of the initial par- 
tons so that we model only the transverse momentum of 
the total ti system, for which we use a prior derived from 
Monte Carlo simulation. This leaves a total of 19 dimen- 
sions over which the integral must be evaluated, which 
we perform using a quasi-Monte Carlo technique. 

Handling of background events is unchanged from our 
previous publication [5]. We identify events likely to 
be background using a jetnet 3.5 artificial neural net- 
work [M] with ten inputs. We construct distributions of 
the neural network output weight u for signal, S(u), and 
background, B(u), events, normalized to their overall ex- 
pected fractions, and calculate the expected background 
fraction for a given event as /bg(u) = B(u)/(B(u) + 
S(u)). 

We calculate the likelihood for all candidate events un- 
der the assumption that they are signal, but the com- 
bined likelihood contains contributions from both signal 
and background events. However, only the signal events 
contain information about m t , so using Monte Carlo- 
simulated events we compute the average likelihood for 
background events and subtract it from the total likeli- 
hood: 

log L adi (m t , A JES ) ( 
= ^2 [log£(j/i|m t , A JES )-f bg (ui) log L bg (m t , A JES )], 

iGevents 

where L a dj is the adjusted total likelihood for a given set 
of events, L(yl\mt, Ajes) is the likelihood for an individ- 
ual event from Eq. ([I]), /bg(uj) is the background fraction 
for a given event with a neural network output Ui, and 
Lbg(nH, Ajes) is the average likelihood for a background 
event. 

Besides background events, the sample includes events 
which contain a real ti, but where one or more of the four 
jets and/or the lepton observed in the detector do not 
come directly from the ti decay, and are not well-modeled 
by the signal likelihood or handled by the background 
subtraction above. These events, which we refer to as 
"bad signal," have a variety of sources (extra jets from 
gluon radiation, ti events where both W bosons decay 
into leptons or hadrons, W — > rv decay, etc.) and make 
up 36% of the simulated ti events for m t — 172.5 GeV/c 2 . 
We suppress these events by requiring that the peak log- 
likelihood value for an event be at least 10. This cut 
retains 96.3% of the signal, while rejecting 30.8% of the 
bad signal and 37.3% of the background. 

We test and calibrate the method by constructing sim- 
ulated experiments using the Monte Carlo samples of ti 
events and background described earlier. For a given in- 
put mt and Ajes, we perform 2000 experiments using a 
Poisson distribution with mean of 1089 events (the num- 
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ber of events expected to pass the likelihood cut), and use 
these to calibrate the measurement as a function of the 
input m t and Ajes- Figure [T] shows the output mass be- 
fore calibration and the calibrated expected uncertainty. 
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FIG. 1: Simulated experiment results using Monte Carlo sig- 
nal and background events. Top: output m t vs. input mt, 
before calibration is applied. Bottom: expected uncertainty 
a m vs. input m t , with calibration applied. The lines are 
linear best fits. 



In the data we find a total of 1087 events which pass 
all of the selection requirements (including the likelihood 
peak cut), of which 854 have 1 b tag and 233 have >1 b 
tag. Figure [2] shows the resulting 2-D likelihood contours 
for 1 cr, 2 a, and 3 a after all calibration. 

CDF Run II 5.6 fb" 1 
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FIG. 2: Measured 2-D likelihood on the data events, with 
the contours corresponding to a 1-a, 2-a, and 3-<r uncertainty 
in the final mt measurement from the profile method. The 
marker shows the point of maximum likelihood. 

To obtain a 1-D likelihood curve in m t only, we treat 
Ajes as a nuisance parameter and eliminate it using the 
profile likelihood method [5T], where we take the maxi- 
mum value of the likelihood along the Ajes axis for each 
m t value. The top quark mass value extracted from the 
profile likelihood after calibration is m t — 173.0 ± 0.9 
GeV/c 2 . We can separate this uncertainty into the sta- 
tistical uncertainty on m t and the uncertainty due to 



Ajes by fixing the Ajes value to its maximum likeli- 
hood value. We find that the uncertainty from the re- 
sulting 1-D likelihood is 0.7 GeV/c 2 , so we assign the 
remaining uncertainty of 0.6 GeV/c 2 to Ajes an d con- 
clude m t = 173.0 ± 0.7 (stat.) ± 0.6 (JES) GeV/c 2 . 

To validate the likelihood cut procedure, we compare 
the peak values of the log-likelihood curves obtained 
with data to those obtained with Monte Carlo-simulated 
events at m,t = 172.5 GeV/c 2 (the nearest available mass 
value) . The results are shown in Fig. [3j 
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FIG. 3: Comparison of the log-likelihood value of the peak 
of likelihood curves for data and Monte Carlo events. The 
vertical line at 10 indicates the likelihood cut used in this 
analysis. A Kolmogorov-Smirnov test gives a confidence level 
of 0.93, showing good agreement between the two. 

The systematic uncertainties on m t , given in Table [TTJ 
are derived using the methods described in Ref. [5], In 
brief, we include uncertainties coming from: the calibra- 
tion method; signal Monte Carlo modeling, evaluated by 
comparing events simulated with the pythia and her- 
WIG [33] generators; variations of the parameters used 
for initial state radiation (ISR) and final state radiation 
(FSR); a residual JES uncertainty because the JES un- 
certainty contains several components with different pt 
and T) dependence; additional uncertainties on the energy 
scale for b jets; uncertainty on the lepton px scale; multi- 
ple hadron interactions, to take into account uncertainty 
on the jet corrections as a function of the number of inter- 
actions in the event; uncertainties arising from the PDFs 
used in the integration; and the background modeling. 
This analysis includes a systematic uncertainty due to 
color reconnection effects, not considered in our previous 
analysis. We use PYTHIA version 6.4.20, which includes a 
color reconnection model [23] , and measure the difference 
between two tunes, Tune A, which is the tune used in this 
analysis, and Tune ACR, which adds color reconnection 
effects to Tune A. The individual systematic uncertain- 
ties are added in quadrature to obtain the final total of 
0.9 GeV/c 2 . 



7 



In conclusion, the measured top quark mass in a sam- 
ple with 5.6 fb _1 of integrated luminosity, with 1087 
events passing all cuts, is m t — 173.0 ± 0.7 (stat.) ± 
0.6 (JES) ± 0.9 (syst.) GeV/c 2 , for a total uncertainty 
of 1.2 GeV/c 2 . The improved integration techniques and 
increased data sample make this the best single measure- 
ment of the top quark mass to date, and it is comparable 
in precision to the most recent combination for the top 
quark mass at the Tevatron [2]. 
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TABLE II: List of systematic uncertainties on mt- 



Systematic source Uncertainty (GeV/c 2 ) 



Calibration 


0.10 


MC generator 


0.37 


ISR and FSR 


0.15 


Residual JES 


0.49 


6- JES 


0.26 


Lepton pt 


0.14 


Multiple hadron interactions 


0.10 


PDFs 


0.14 


Background modeling 


0.33 


Color reconnection 


0.37 


Total 


0.88 



